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Astrophysical observations spanning dwarf galaxies to galaxy clusters indicate that dark matter (DM) halos 
are less dense in their central regions compared to expectations from collisionless DM N-body simulations. 
Using detailed fits to DM halos of galaxies and clusters, we show that self-interacting DM (SIDM) may provide 
a consistent solution to the DM deficit problem across all scales, even though individual systems exhibit a wide 
diversity in halo properties. Since the characteristic velocity of DM particles varies across these systems, we 
are able to measure the self-interaction cross section as a function of kinetic energy and thereby deduce the 
SIDM particle physics model parameters. Our results prefer a mildly velocity-dependent cross section, from 
o/m ~ 2 cm 2 /g on galaxy scales to a/m ~ 0.1 cm 2 /g on cluster scales, consistent with the upper limits 
from merging clusters. Our results dramatically improve the constraints on SIDM models and may allow the 
masses of both DM and dark mediator particles to be measured even if the dark sector is completely hidden 
from the Standard Model, which we illustrate for the dark photon model. 


I. Introduction. Dark matter (DM) is the dominant form 
of matter in galaxies and galaxy clusters, influencing the ob¬ 
served motions of stars and gas, as well as the lensing of light 
rays from distant sources. In turn, these observables provide 
gravitational tracers of DM’s mass distribution, which can 
shed light on its unknown microphysical properties. One im¬ 
portant question is whether DM particles undergo elastic scat¬ 
tering interactions with one another. This process, known as 
DM self-interactions, can leave an imprint on the stmcture of 
DM halos ID. Self-interactions lead to heat transport between 
the inner and outer parts of the halo, typically causing the in¬ 
ner halo to become less dense. Numerical simulations have 
shown that self-interacting DM (SIDM) halos have reduced 
central densities compared to the cold DM (CDM) paradigm, 
where DM is collisionless and non-interacting EH3- 

There are long-standing anomalies from observations of 
dwarf galaxies that their inner regions have less DM mass 
compared to CDM predictions |[8][9]. More recent astrophys¬ 
ical measurements have shown this DM mass deficit problem 
to be ubiquitous: it appears in satellite galaxies of the Milky 
Way and Andromeda HMD, field dwarf galaxies within the 
Local Group EKH, low surface brightness (LSB) spiral 
galaxies HME], high surface brightness spiral galaxies d). 
and galaxy clusters ca. Different astrophysical mechanisms 
have been proposed to address these issues, including feed¬ 
back from supernovae l20)[2Tll and active galactic nuclei l22l . 
However, each mechanism cannot by itself solve the mass 
deficit problem everywhere. 

The SIDM paradigm has the potential to solve the mass 
deficit problem in halos across all scales. However, a quan¬ 
titative comparison with observational data is lacking, in par¬ 
ticular for spiral galaxies and clusters. Previous studies relied 
on phenomenological profiles that share qualitatively similar 
features as SIDM — namely, “cored” profiles with reduced 
central densities compared to CDM — but it is unclear how 


closely these mimic true SIDM profiles inferred from sim¬ 
ulations. It is also unknown what inferences to draw about 
the self-interaction cross section and whether SIDM can suc¬ 
cessfully accommodate the diversity of DM halos exhibited in 
nature, with core sizes from dwarf to cluster scales ranging 
~ 0.5 — 50 kpc. Moreover, cored profiles do not capture the 
gravitational effect from baryons in regions that are baryon- 
dominated, as in the central regions of clusters, which tends 
to increase the DM density l23l . 

In this Letter , we address these important issues. In Sec. II, 
we present a semi-analytical method for modeling SIDM ha¬ 
los S3 [23]. While numerical simulations remain the standard 
for studies of DM structure, it is not yet feasible for them to 
map out SIDM parameter space because even the simplest 
particle models have a rich diversity of velocity-dependent 
cross sections Il24l l25il. Our halo model provides an important 
bridge between simulations, astrophysical data, and particle 
theories for SIDM. In Sec. Ill, we fit our model to a sample 
of twelve dwarf and LSB galaxies and six clusters of galaxies 
and infer the self-interaction cross section for these systems. 
We also test our method against mock data from N-body sim¬ 
ulations of SIDM model to show that we are able to reproduce 
the correct input cross section. In Sec. IV, we discuss the role 
of halo concentration in explaining the observed diversity of 
halo profiles for DM-dominated galaxies. 

The key insight resulting from our analysis is that data on 
different scales, from dwarf galaxies to clusters, provide com¬ 
plementary handles to disentangle DM microphysics. Since 
DM velocity dispersion increases with halo mass, studying ha¬ 
los from dwarf galaxies to galaxy clusters allows one to probe 
the self-interaction cross section as a function of kinetic en- 
ergy — similar to tuning the beam energy in a particle collider. 
In Sec. V, we illustrate this key point by applying our results 
to a minimal dark photon model where self-interactions are 
described by a Yukawa potential. In Sec. VI, we summarize 
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our main conclusions. 

II. SIDM halo model. Scattering between DM particles 
is more prevalent in the halo center where the DM density is 
largest. It is useful to divide the halo into two regions, sepa¬ 
rated by a characteristic radius rq where the average scatter¬ 
ing rate per particle times the halo age (i age ) is equal to unity. 
Thus, 


rate x time 
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where a is the scattering cross section, m is the DM parti¬ 
cle mass, v is the relative velocity between DM particles and 
(...) denotes ensemble averaging. Since we do not assume 
cr to be constant in velocity, we find it more convenient to 
quote ( av) /m rather than a /to. We set i age = 5 and 10 Gyr 
for clusters and galaxies, respectively. Although Eq. 0 is a 
dramatic simplification for time integration over the assembly 
history of a halo, we show by comparing to numerical simu¬ 
lations that it works remarkably well. 

For halo radius r > r i, where scattering has occurred 
less than once per particle on average, we expect the DM 
density to be close to a Navarro-Frenk-White (NFW) profile 
p(r) = p s (r/r s )~ 1 (l+r/r s )~ 2 characteristic of collisionless 
CDM l26l . In the halo center, for radius r < r i, scattering 
has occurred more than once per particle. Here, we expect 
DM particles to behave like an isothermal gas satisfying the 
ideal gas law p = pa q, where p , p are the DM pressure and 
mass density and ao is the one-dimensional velocity disper¬ 
sion. Since the inner halo achieves kinetic equilibrium due 
to DM self-interactions, the density profile can be determined 
by requiring hydrostatic equilibrium, Vp = — pV'ktot- Here, 
<f>tot is the total gravitational potential from DM and bary- 
onic matter, which satisfies Poisson’s equation V 2< f>tot = 
4-7 tG(p + pb), where G is Newton’s constant and pb is the 
baryonic mass density. These equations yield 

(7q V 2 In p = -47rG(p + Pb) , (2) 


which we solve to obtain p(r) assuming spherical symmetry. 

We model the full SIDM profile by joining the isothermal 
and collisionless NFW profiles together at r = r\\ 


( Piso(r) , r < rq 
\ PNFw(r) , r > r i 
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where pi so is the solution to Eq. 0- We fix the NFW param¬ 
eters ( p s , r s ) by requiring that the DM density and enclosed 
mass for the isothermal and NFW profiles match at rq. Thus, 
our SIDM halo profile is specified by three parameters: the 
central DM density p 0 = p(0), velocity dispersion a 0 , and 
rq. Lastly, we note that this model exhibits a two-fold degen¬ 
eracy in solutions for (av)/m. We keep the smaller (av)/m 
solutions but note that this situation may be indicative of the 
degeneracy between halo profiles with cores that are growing 
or shrinking in time 0. 

III. SIDM fits. To constrain DM self-interactions, we con¬ 
sider a set of six relaxed clusters and twelve galaxies with 



FIG. 1: Self-interaction cross section measured from astrophysical 
data, given as the velocity-weighted cross section per unit mass as 
a function of mean collision velocity. Data includes dwarfs (red), 
LSBs (blue) and clusters (green), as well as halos from SIDM 
N-body simulations with cr/m = 1 cm 2 /g (gray). Diagonal 
lines are contours of constant a/m and the dashed curve is the 
velocity-dependent cross section from our best-fit dark photon model 
(Sec. V). 


halo masses spanning 10 9 — 10 15 M 0 . These objects ex¬ 
hibit central density profiles that are systematically shallower 
than p oc r _1 predicted from CDM simulations. To determine 
the DM profile for each system, we perform a Markov Chain 
Monte Carlo (MCMC) scan over the parameters (p 0 ,cro,rq) 
characterizing the SIDM halo, as well as the mass-to-light ra¬ 
tio Y, for the stellar density. The value for p(j \) determines 
the velocity-weighted cross section (av) /m from Eq. 0, as a 
function of average collision velocity ( v) = (4/y/ir)ao for 
a Maxwellian distribution. We also verify our model and 
MCMC fit procedure using a mock data set from simulations. 

Clusters. We consider the relaxed clusters from the data 
set of Newman, et al. lfl9l [27] for which spherical modeling 
is appropriate (MS2137, A611, A963, A2537, A2667, and 
A2390). These clusters have stellar kinematics as well as 
strong and weak lensing measurements allowing the mass pro¬ 
file to be measured from stellar-dominated inner region (~ 10 
kpc) out to the virial radius (~ 3 Mpc). The baryonic and 
DM densities are disentangled by constraining T* through 
the assumption that all the clusters share a similar star for¬ 
mation history. The inferred DM density profile is consistent 
with CDM expectations except in the inner Of 10) kpc region 
where a mass deficit is inferred ED- These small core sizes 
dictate the preference for a velocity-dependent cross section. 

We model each cluster using Eq. 0 and fit directly to the 
stellar line-of-sight velocity dispersion data 127). We include 
the gravitational effect of the stars following Eq. 0 and allow 
for a ±0.1 dex spread in pb to account for systematic uncer¬ 
tainties HU [271. Further, as a proxy for fitting to the gravi¬ 
tational lensing data at large radii, we fit to posteriors of the 
maximum circular velocity T4 max and the corresponding radius 
v max that have been obtained from the lensing data ED. 
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FIG. 2: Top: SIDM density profile fit to cluster A2537 (orange) com¬ 
pared to NFW profile (cyan) and comparison to stellar kinematics 
data (inset). Bottom: SIDM fit to the rotation curve of galaxy IC2574 
(orange) with contributions from the SIDM halo (solid), the gas disk 
(dashed), and stellar disk (dotted). 


As an example, we show our results for cluster A2537 in 
Fig. |2](7op). Our SIDM lit is shown by the orange band (la- 
width) and the dashed line shows the mean. The CDM predic¬ 
tion (cyan) is the NFW profile obtained from the gravitational 
lensing data G3, which provides a poor fit to the stellar kine¬ 
matic data (red boxes in inset figure). The black point is the 
value of n and its lo- width. It is reassuring that the CDM and 
SIDM fits, while agreeing at large radii, begin to diverge at n. 
The inferred values of ( av)/m for all six clusters are shown 
in Fig. [T] (green points). Fitted with a constant cross section, 
we find a/m = 0 . 10 ^o.o 2 cm 2 /g. 

Dwarf and Low Surface Brightness Galaxies. To mea¬ 
sure DM self-interactions at small-to-intermediate scales, we 
consider rotation curves of five dwarf galaxies (IC 2574, NGC 
2366, Ho II, M81 dwB, DDO 154) in the THINGS sam¬ 
ple d and seven LSB galaxies (UGC 4325, F563-V2, F563- 
1, F568-3, UGC 5750, F583-4, F583-1) from Kuzio de Naray, 
et al. (29l . Two galaxies have been omitted from each of these 
samples for which I4iax was not well-determined. 

To model these galaxies, we include the contributions to 
the rotation curve from DM, gas, and stars, with T, allowed 
to vary uniformly by ±0.3 dex from the quoted population 


synthesis values |[28l [30l l. We have checked that it is a good 
approximation to neglect the gravitational effect of baryons on 
the SIDM density profile in Eq. |2|. In our likelihood, we also 
include a systematic error (in quadrature with the statistical 
error) of 5% of the last measured velocity to avoid skewing 
our fits based on data points with small errors, (7(lkm/s), 
since non-circular motions cannot be excluded at this level. 

As an example, we show the SIDM fit to the rotation curve 
of IC2574 in Fig. ^(Bottom). The inferred values of (av)/m 
for the galaxies, shown in Fig. [I] evidently prefer a larger a/m 
than the cluster measurement. Fitting all twelve galaxies with 
a constant cross section, we find a/m = 1.9^4 cm 2 /g. We 
note that this value does not include systematic errors, which 
we discuss next. 

Simulated halos. To test our analytic model, we created 
mock rotation curve data from halos in a/m = 1 cm 2 /g sim¬ 
ulations (without baryons) and fit them with our model. Each 
rotation curve consisted of 20 points with a uniform 10% ve¬ 
locity error and covering a range 0.1 < r/r s < 3 . We chose 
six halos with virial masses in the range 10 11 —10 14 Mq from 
Ref. 0 and two dwarf-sized halos around 10 10 AT., from 
Ref. Q. 

The fit results shown by the gray points in Fig. |T] demon¬ 
strate that our simple halo model is in good agreement with 
results from cosmological N-body simulations for SIDM, ex¬ 
cept for the presence of a bias toward larger cross sections by 
a factor of ~ 2. The open circles, which also line up along 
a/m = 1 cm 2 /g, represent our SIDM profiles matched onto 
the “true” NFW profile for the same halos simulated without 
DM self-interactions 00. This analysis supports the sim¬ 
ple picture in our model that the SIDM halo properties may 
be approximated by the corresponding CDM halo properties 
augmented with a core determined by Eq. Q. 

IV. Diversity. There is considerable diversity in the prop¬ 
erties of the galaxy cores, with almost an order of magnitude 
spread in density at fixed U max ED. This has also been re¬ 
cently emphasized in terms of 14(2 kpc), the measured circu¬ 
lar velocity at 2 kpc ED , which shows a factor of 2 — 3 scatter 
for halos with 50km/s < V max < lOOkm/s. This diversity 
is also reflected in the scatter in central values for (av) /m for 
the galaxies in Fig. |T] 

How does this scatter arise in our model? The answer is 
surprising in its simplicity: it is directly related to the halo 
assembly history. Different formation histories encoded in 
( p s , r s ) values (essentially the CDM halo-to-halo scatter) lead 
to SIDM halos with different core sizes and central densities 
through Eq. 0. This explanation is implicit in Fig. [I] where 
the large errors on (av) reflect, partly, the lack of constraints 
on ( p s ,r s ). Choosing the “right” value of ( p s ,r s ) for each 
galaxy would reduce the scatter in (av)/m considerably. 

If we fix the p s -r s relation to its median in ACDM cos¬ 
mology f32l in our analysis, the galaxies UGC 5750 and IC 
2574 prefer the largest cross sections, a/m ~ 10 cm 2 /g, 
while M81 dwB prefers the smallest cross sections, a/m ~ 
0.1 cm 2 /g. However, if UGC 5750 and IC 2574 halos are 2a 
less concentrated and M81 dwB halo 2 a more concentrated 
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than the median halo, the inferred cr/m become consistent, 
within errors, with ~ 1 cm 2 /g. 

Turning this around, we can fix <t/to and look at the impact 
of the scatter in the p s -r s relation on 14(2 kpc). Within our 
analytic model, we have checked that the spread in the p s - 
r s relation in ACDM leads to about a factor of two spread in 
14(2 kpc) for the relevant galaxies. If we were to add baryons 
(which could be important within 2 kpc), it is conceivable that 
the bulk of the spread seen in Ref. PTl can be explained. 

V. Dark matter particle properties. The energy depen¬ 
dence of the cross section allows one to discern the underly¬ 
ing particle dynamics of SIDM. The data in Fig. [[[range over 
a factor of 10 4 in kinetic energy and prefer a cross section that 
mildly falls with energy. 

To illustrate the implications for particle physics, let us con¬ 
sider the dark photon model for DM self-interactions. In this 
model, self-interactions are governed by a Yukawa potential, 
V(r) = a'e~^ r /r, where a' is the coupling constant (anal¬ 
ogous to the fine structure constant a « 1/137) and p is the 
dark photon mass, which screens the potential 13314351 . To be 
concrete, we will set a' = a. We then compute (av)/m using 
standard partial wave methods discussed in Ref. l25l . Com¬ 
paring the theoretical predictions to the data points in Fig. [T] 
using a Ay 2 test, we determine the preferred regions for the 
DM mass m and dark photon mass p . To take into account 
the uncertainty in our modeling (apparent in our predictions 
for the simulated halos), we have included an additional sys¬ 
tematic uncertainty (in quadrature) of A (log (crt;)/m) = 0.3 
and A(log(w)) = 0.1 for each system. 

Our results shown in Fig. [3]illustrate the important comple¬ 
mentarity between observations across different scales in con¬ 
straining DM microphysics. The red, blue, and green shaded 
bands show the individual 95% confidence level (CL) regions 
preferred by our analysis of dwarf galaxies, LSBs, and clus¬ 
ters, respectively. The solid (dashed) black contour shows the 
95% (99%) CL region from all observations combined. These 
data prefer DM mass of 15^4 GeV and dark photon mass of 
17 ± 4 MeV at 95% CL. For the best-fit values of m and p, 
we plot (av) /m as a function of (v) in Fig.[l](dashed). 

Fig. [3] also shows the regions excluded by the Bullet Clus¬ 
ter constraint of a/m < 1.25 cm 2 /g at 68% CL |[36l at 
v = 2000 km/s (dot-dashed) and the constraint from an en¬ 
semble of merging clusters of a/m < 0.47 cm 2 /g at 95% 
CL ll37l at v = 900 km/s (long-dashed). A more refined anal¬ 
ysis of the merging clusters, including large dissociative clus¬ 
ters that show offsets between the luminous and dark compo¬ 
nents lf36l[384 - 14011 . would be interesting in light of the velocity 
dependence. 

It is remarkable that astrophysical observations can pick out 
a closed range for the DM mass m (albeit within the sim¬ 
ple model we have adopted). For to < 10 GeV, the cross 
section changes little with velocity, a oc m 2 / p 4 , which is 
disfavored by the velocity dependence evident in Fig. [T] For 
to )Z 100 GeV, cross section tends to the Rutherford limit, 
a oc 1/(m 2 v 4 ), which is too steep a velocity dependence to be 
consistent with our fits. The preferred region lies in between 
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FIG. 3: Parameter space for dark photon model of self-interactions 
(with a' = a), preferred by dwarfs (red), LSB spiral galaxies (blue), 
and clusters (green), each at 95% CL. Combined 95% (99%) region 
is shown by the solid (dashed) contours. The estimated Bullet Clus¬ 
ter excluded region lies below dot-dashed curve and the ensemble 
merging cluster excluded region below the long-dashed curve. 


these extremes: a is constant at small velocity and turns over 
to a Rutherford-like dependence at large velocity. 

VI. Conclusions. SIDM paradigm may provide a unified 
explanation for the apparent deficit of DM in the central re¬ 
gions of galaxies and clusters. We have explored the direct 
connection between self-interactions and astrophysical obser¬ 
vations for a set of twelve galaxies and six clusters using a 
simple model for SIDM halos calibrated to N-body simula¬ 
tions. Despite the diversity of DM halo properties in these 
systems, the majority of dwarfs and LSBs is remarkably con¬ 
sistent with a/m « 2 cm 2 /s. Clusters favor 0.1 cm 2 /g 
because their halo profiles are largely consistent with CDM 
except in the inner 0(10kpc) region. The velocity depen¬ 
dence discernible in these data provides an important step to¬ 
ward understanding the possible particle physics of DM self¬ 
interactions. Within the dark photon model we considered, 
these data prefer DM mass of ~ 15 GeV and dark photon 
mass of ~ 17 MeV. While these conclusions are model- 
specific, SIDM in general indicates a new mass scale much 
below than the electroweak scale. Using DM halos as particle 
colliders, we may be able to unveil the particle physics nature 
of DM, independent of whether the dark and visible sectors 
are coupled via interactions beyond gravity. 
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